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Abstract 


The present article discusses some interesting phenomena including the Lense-Thirring type anomalous 
precession, using a known spherical kinetic dynamics approach. Other implications include a plausible 
revised version of the celestial quantization equation described by Nottale and Rubcic & Rubcic. If the 
proposition described herein corresponds to the facts, then this kinetic dynamics interpretation of 
‘frame-dragging’ effect could be viewed as a step to unification between GTR-type phenomena and 
QM. Further observation to verify or refute this conjecture is recommended, plausibly using LAGEOS- 
type satellites. 
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Introduction 


It is known, that the use of Bohr radius formula to predict celestial quantization has led to numerous verified 
observations [1]. This approach was based on Bohr-Sommerfeld quantization rules [2][3]. While this kind of 
approach is not widely accepted yet, this could be related to wave mechanics equation to describe large-scale 
structure of the Universe [4], and also a recent suggestion to reconsider Sommerfeld’s conjectures in Quantum 
Mechanics [5]. Some implications of this quantum-like approach include exoplanet prediction, which becomes a 
rapidly developing subject in recent years [6][7]. 

Rubcic & Rubcic’s approach [2] is particularly interesting in this regard, because they begin with a conjecture 
that Planck mass (m, = Vvic/2zG ) is the basic entity of Nature, which apparently corresponds to Winterberg’s 
assertion of superfluid Planckian aether comprised of phonon-roton pairs [8]. In each of these pairs, superfluid 


vortices can form with circulation quantized according to fv, .dx =nh/m,. This condition implies the 


Helmholtz vortex theorem, d / athy,. .dx = 0 . This relationship seems conceivable, at least from the viewpoint 


of likely neat linkage between cosmology phenomena and various low-temperature condensed matter physics 
[9][10][11]. In effect, celestial objects at various scales could be regarded as spinning Bose-Einstein condensate; 
which method has been used for neutron stars [32]. 

Despite these aforementioned advantages, it is also known that all of the existing celestial quantization 
methods [1][2][3] thus far have similarity that they assume a circular motion, while the actual celestial orbits (and 
also molecular orbits) are elliptical. Historically, this was the basis of Sommerfeld’s argument in contrast to 
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Bohr’s model, which also first suggested that any excess gravitational-type force would induce a precessed orbit. 
This is the starting premise of the present article, albeit for brevity we will not introduce elliptical effect yet [12]. 

Using a known spherical kinetic dynamics approach, some interesting phenomena are explained, including 
the receding Moon, the receding Earth from the Sun, and also anomalous precession of the first planet (Lense- 
Thirring effect). Despite some recent attempts to rule out the gravitational quadrupole moment (J2) contribution 
to this effect [13][14][15][16][17], it seems that the role of spherical kinetic dynamics [12] to Lense-Thirring 
effect has not been taken into consideration thus far, at least to this author’s knowledge. 

After deriving prediction for some known observed phenomena, this article will also present a revised version 
of quantization equation of L. Nottale [1] in order to take into consideration this spherical kinetic dynamics effect. 
If the proposition as described here corresponds to the facts, then this approach could be viewed as a step to 
unification of GTR-type phenomena and Quantum Mechanics. 


Spherical kinetic dynamics, Earth bulging effect 


In this section we start with some basic equations that will be used throughout the present article. It is 
assumed that the solar nebula is disk-shaped and is in hydrostatic equilibrium in the vertical direction. Let 
suppose that the disk has approximately Keplerian rotation, œ; then the half-thickness of the disk is given by 

[4d; p.4-5]: 


d=c,/a@ (1) 
and 


c, ¥vVkT/m (la) 


where d and cs represents half-thickness of the disk and sound velocity, respectively. 
In order to find the spherical kinetic dynamics contribution to Lense-Thirring effect, we begin with the 
spinning dynamics of solid sphere with mass M. Using the known expression [12; p.6, p.8]: 


Epines = -10/2 (2) 


kinetic 
I =2MR°/5 (2a) 


sphere 


where Iz, œ, M, R represents angular momentum, angular velocity, spinning mass of the spherical body, and 
radius of the spherical body, respectively. Inserting equation (2a) into (2) yields: 


E....=—-MR’*o@* /5 (3) 


kinetic 
This known equation is normally interpreted as the amount of energy required by a spherical body to do its 
axial rotation. But if instead we conjecture that ‘galaxies get their angular momentum from the global rotation of 
the Universe due to the conservation of the angular momentum’ [34], and likewise the solar system rotates 
because of the corresponding galaxy rotates, then this equation implies that the rotation itself exhibits extra kinetic 
energy. Furthermore, it has been argued that the global rotation gives a natural explanation of the empirical 


relation between the angular momentum and mass of galaxies: J ~ aM °”? [34]. This conjecture is also relevant 
in the context of Cartan torsion description of the Universe [18]. For reference purpose, it is worthnoting in this 
regard that sometime ago R. Forward has used an argument of non-Newtonian gravitation force of this kind, 
though in the framework of GTR (Amer.J.Phys. 31 No. 3, 166, 1963). 

Let suppose this kind of extra kinetic energy could be transformed into mass using a known expression in 
condensed-matter physics [10b; p.4], with exception that cs is used here instead of v to represent the sound 
velocity: 


E inetic Ms p) = C,.p = m, e (4) 
where the sound velocity obeying [10b; p.4]: 
c,’(n) =(n/m)(d? €/ dn’) (4a) 


Physical mechanism of this kind of mass-energy transformation is beyond the scope of the present article, 
albeit there are some recent articles suggesting that such a condensed-matter radiation is permitted [35]. Now 
inserting this equation (4) into (3), and by dividing both sides of equation (3) by Ar, then we get the incremental 
mass-energy equivalent relation of the spinning mass: 


Am, / At = —.(A@/ At).MR? [(5.c,”) (5) 
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By denoting © = Aw@/ At , then this equation (5) can be rewritten as: 
AM / At = -@.MR’al(5.c,’) (5a) 
For © = 0 the equation (Sa) shall equal to zero, therefore this equation (Sa) essentially says that a linear change 
of angular velocity observed at the surface of the spinning mass corresponds to mass flux, albeit this effect is 
almost negligible in daily experience. But for celestial mechanics, this effect could be measurable. 


If, for instance, we use the observed anomalous decceleration rate [30] of angular velocity of the Earth as 
noted by Kip Thorne [19]: 


lo /o| = 6x10" years (6) 
And using values as described in Table 1 for other parameters: 


Table 1. Parameter values to compute kinetic expansion of the Earth 























Parameter Value Unit 
Re 6.38x10° m 
Me 5.98x1074 kg 
Te 2.07x10° sec 
We 3.04x10° rad/s 
Cs 0.14112 m/s 














It is perhaps worth noting that the only free parameter here is cs =0.14112 m/sec. This value is approximately 
within the range of Barcelo et al.’s estimate of sound velocity (at the order of cm/sec) for gravitational Bose- 
Einstein condensate [11], provided the Earth could be regarded as a spinning Bose-Einstein condensate. 
Alternatively, the sound velocity could be calculated using equation (la), but this obviously introduces another 
kind of uncertainty in the form of determining temperature (T) inside the center of the Earth; therefore this method 
is not used here. 

Then inserting these values from equation (6) and Table 1 into equation (5a) yields: 

AM / At = 3.76x10'°kg/ year (7) 

Perhaps this effect could be related to a recent Earth bulging data, which phenomenon lacks a coherent 
explanation thus far [36]. 


Prediction of the receding Moon and the receding planets from the Sun 


Now let suppose this predicted value (7) is fully conserved to become inertial mass, and then we could 
rewrite Nottale’s method of celestial quantization [1]. Alternatively, we could begin with the known Bohr- 
Sommerfeld quantization rule [3]: 


f p,da, = n,.2m.e I(a,.c) (8a) 
Then, supposing that the following substitution is plausible [3]: 
ela, >GMm/a, (8b) 


where e,Qe,2 represents electron charge, Sommerfeld’s fine structure constant, and gravitational-analogue of fine 
structure constant, respectively. This corresponds to Nottale’s basic equations v, =a@,.c/n=v,/n and vo=144 
km/sec [1]. And by introducing the gravitational potential energy [12]: 
®(r,.9)=-GM 1 r1- J,.(a1r} Bco? 8—1)/2| (8c) 
where Qis the polar angle (collatude) in spherical coordinate, M the total mass, and a the equatorial radius of 
the solid. 
Neglecting higher order effects of the gravitational quadrupole moment J2 [13][14][15][16][17], then we get 
the known Newtonian gravitational potential: 
@®=-GM/r (8d) 
Then it follows that the semi-major axes of the celestial orbits are given by [1][3]: 
r, = GMn* lv, (8e) where n=1,2,....is the principal quantum 
number. 
It could be shown, that equation (8a) also corresponds to the conjecture of quantization of circulation [4b]. 
By reexpressing equation (8e) for mass flux effect (5) by defining M,,,=M,+AM,/At,, then the total 
equation of motion becomes: 
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(M +AM /At) =(r+Ar/At).v) (Gn?) (8f) 
For A —> 0, equation (8f) can be rewritten as: 


dM / dt — y.dr/dt+M -r.y =0 (8g) 
where 
=v Gn’) (8h) 
Now inserting (5a) into equation (8g), and dividing both sides by y, yields: 
dr/dt-M/y+r+@MR’ol(75.c, ) =0 (8i) 
This equation (81) can be rewritten in the form: 
r+r+@o=0 (8j) 
by denoting ř =dr/dt and 
p=-M/y[1-aR’ol5.c,)] (8k) 


if we suppose a linear decceleration at the surface of the spinning mass. Equation (8j) and (8k) is obviously a first- 
order linear ODE equation [26], which admits exponential solution. In effect, this implies that the revised equation 
for celestial quantization [1][2] takes the form of spiral motion. This could also be interpreted as a plausible 
solution of diffusion equation in dissipative medium [33], which perhaps may also correspond to the origin of 
spiral galaxies formation [28]. And if this corresponds to the fact, then it could be expected that the spiral galaxies 
and other gravitational clustering phenomena [22b] could also be modeled using the same quantization method 
[39], as described by Nottale [1] and Rubcic & Rubcic [2]. 

To this author’s knowledge these equations (8J) and (8k) have not been presented before elsewhere, at least in 
the context of celestial quantization. 

Inserting result in equation (7) into (8e) by using n=3 and vo=23.71 km/sec for the Moon [2] yields a receding 
orbit radius of the Moon as large as 0.0401 m/year, which is very near to the observed value ~ 0.04 m/year [20]. 
The quantum number and specific velocity here are also free parameters, but they have less effect because these 
could be replaced by the actual Moon orbital velocity using v, =v, /n [1]. 


While this kind of receding Moon observation could be described alternatively using oscillation of 
gravitational potential [30], it seems the kinetic expansion explanation is more preferable particularly with regard 
to a known phenomenon of continental drift [29], and also perhaps a known periodic geological layering described 
by Alvarez et al. Furthermore, using this kind of approach it seems that we could also offer a plausible explanation 
for the kinetic origin of volcanoes eruption. Apparently, none of these effects could be explained using oscillation 
of gravitational field argument, because they are relentless effects. 

In this regard, it is interesting to note that Sidharth has argued in favor of varying G [21]. From this starting 
point, he was able to explain —among other things-- anomalous precession (Lense-Thirring effect) of the first 
planet and also anomalous Pioneer acceleration. This will be discussed in the subsequent section. In principle, 
Sidharth’s basic assertion is [21]: 

G=G,{1+t/ts) (9) 

It is worthnoting here that Barrow [40c] has also considered a somewhat similar argument in the context of 
varying constants: 

G=Gytg/(t—c) (9a) 

However, in this article we will use (9) instead of (9a), partly because it will lead to more consistent predictions 
with observation data. Alternatively, we could also hypothesize using Maclaurin formula: 

G=G,.e'? =G- (l+t/to t(t/te) /24(t/tey 134...) (9b) 

This expression is a bit more consistent with the exponential solution of equation (8J) and (8k). Therefore, 
from this viewpoint equation (9) could be viewed as first-order approximation of (9b), by neglecting second and 
higher orders in the series. It will be shown in subsequent sections, that equation (9) is more convenient for 
deriving predictions. 

If we conjecture that instead of varying G, the spinning mass M varies, then it would result in the same effect 
as explained by Sidharth [21], because for Keplerian dynamics we could assert k=GM, where k represents the 
stiffness coefficient of the system. Accordingly, Gibson [22] has derived similar conjecture of exponential mass 
flux from Navier-Stokes gravitational equation, which can be rewritten in the form: 

M =G,.e'® =M,.+t/tyt(t/te) /24(t/te) 134...) (9c) 
provided we denote for consistency [22]: 
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to =T, 27 (9d) 

Using the above argument of Maclaurin series, equation (9c) could be rewritten in the similar form with 
(9) by neglecting higher order effects: 

M =Mo(l+t/tə) (10) 

In a recent article Gibson & Schild [23] argue that their gravitational Navier-Stokes approach results in 
better explanation than what is offered by Jeans instability. Furthermore, R.M. Kiehn has also shown that the 
Navier-Stokes equation corresponds exactly to Schroedinger equation [27], which seems to support the idea 
of quantization of celestial motion [1][2][3]. A plausible extension of Euler equation and Jeans instability to 
describe gravitational clustering has been discussed in [22b], which corresponds to viscosity term and also 
turbulence phenomena [22c,22d] described by Gibson. Therefore, apparently equation (10) is more 
consistent with kinematical gravitational instability consideration than (9). 

From equation (10) we could write for M at time difference At =t, -t : 


M, =M,A1+t,/ta) (11) 

M, =M,Al+t,/ta) (12) 
from which we get: 

AM =(M ,/te)-(t, —t,) (13) 
Inserting our definition At =t, —t, yields: 

AM / At = (M 9/tg) =k (14) 


For verification of this assertion, we could use equation (14) instead of (Sa) to predict mass flux of the Earth. 
Inserting the present mass of the Earth from Table 1 and a known estimate of Earth epoch of 2.2x10° years, we 
get k=0.272x10"° kg/year, which is approximately at the same order of magnitude (ratio=13.83) with equation 
(7). 


Inserting equation (14) into equation (Sa), we get: 


Molts ~—@.MR*al(5.c’) (15) 

which is the basic conjecture of the present article. From this viewpoint we could rewrite equation (8j) and (8k): 
dr /dt =M/y[l-o.R’ol(5.c,)\—r (15a) 

and inserting equation (15), we get: 
dr/dt=M/y..1+M,MtgM)\|—-r (15b) 


A plausible test of this conjecture could be made by inserting this result (14) into equation (8e) and using 
Mo =1.98951x10” g and tẹ = 2x10" year as the epoch of the solar system [21], and specific velocity v.=144 


km/sec [1], then from equation (15b) we get a receding orbit radius for Earth at the order of: 
Ar,,,,/ At = 6.03m/ year (16) 


Earth 

Interestingly, there is an article [24] hypothesizing that there is a tad effect of receding Earth orbit from the 
Sun at the order of 7.5 m/year, supposing Earth orbit radius has been expanding as large as 93x10° miles since 
the beginning of the solar epoch at tẹ = 2x10'° year ago (in the quoted article, it was assumed that the epoch is 
4.5x10° years). Of course, it shall be noted that there is large uncertainty of the estimate of solar epoch, for instance 
Gibson prefers 4.6x10"” sec (or 1.46x10!°year, see [22]). Therefore, it is suggested here to verify this assumption 
of solar epoch using the same tad effect for other planets. For observation purposes, some estimate values were 
presented in Table 2 using the same approach with (15b). 


Table 2. Prediction of planetary orbit radii (r) increment 

















Celestial object Quantum number (n) Orbit increment (m/yr) 
Mercury 3 2.17 
Venus 4 3.86 
Earth 5 6.03 
Mars 6 8.68 




















Quantization of anomalous celestial precession 


It is known that the Newtonian gravitation potential equation (8d) is only weak-field approximation, and 
that GTR makes a basic assertion that this equation is exact. And if the gravitation could be related to boson 
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condensation phenomena [9][10][11], then it seems worth to quote a remark by Consoli [9b; p.2]: “for weak 
gravitational fields, the classical tests of general relativity would be fulfilled in any theory that incorporates 
the Equivalence Principle.’ And also [9b; p.18]: “Einstein had to start from the peculiar properties of 
Newtonian gravity to get the basic idea of transforming the classical effects of this type of interaction into a 
metric structure. For this reason, classical general relativity cannot be considered a dynamical explanation of 
the origin of gravitational forces.” Furthermore, Consoli also argued that the classical GTR effects other than 
anomalous precession could be explained without introducing non-flat metric, as described by Schiff [9b; 
p.19], therefore it seems that the only remarkable observational ‘proof’ of GTR is anomalous precession of 
the first planet [37]. Therefore, it seems reasonable to expect that the anomalous precession effect could be 
predicted without invoking non-flat metric, which suggestion is particularly attributed to R. Feynman, who 
‘believed that the geometric interpretation of gravity beyond what is necessary for special relativity is not 
essential in physics.’ [9d] It will be shown that a consistent approach with equation (10) will yield not only 
the anomalous celestial precession, but also a conjecture that such an anomalous precession is quantized. 
By using the same method as described by Sidharth [21], except that we assert varying mass M instead of 
varying G — in accordance with Gibson’s solution [22]--, and denoting the average angular velocity of the planet 


by 

Q=22/T (17) 
and period T, according to Kepler’s Third Law: 

T =22.a*” /VGM (18) 
Then from equation (10), (17), (18) we get: 

Q-Q, =-Q, t/t, (19) 
Integrating equation (19) yields: 

o(t)=Q-Q, =-(a!1T)t It (20) 


which is average precession at time ‘t’. Therefore the anomalous precession corresponds to the epoch of the 
corresponding system. For Mercury, with T=0.25 year, equation (20) yields the average precession per year 
at time ‘t’: 


TD) merus =Q- Oy = 4rd’ It, (21) 


Using again t = 2x10" year as the epoch of the solar system and integrating for years n=1 ... 100, equation 


(21) will result in total anomalous precession in a century: 
n=100 
a(n) = Y a(n) = 43.86"' percentury 


n=l 


(22) 


It would be more interesting in this regard if we also get prediction of this effect for other planets using the 
same method (20), and then compare the results with GTR-Lense-Thirring prediction. Table 3 presents the result, 


in contrast with observation by Hall and also prediction by Newcomb, which are supposed to be the same [25]. 
Table 3. Comparison of prediction and observed anomalous precession 
































Celestial | Period, O. Hall/ Diff. GTR/ Diff. 
Object T prediction | Newcomb Thirring 

(year) (arcsec/cy) | (arcsec/cy) (%) (arcsec/cy) (%) 
Mercury | 0.25 43.86 43.00 2.03 42.99 -0.05 
Venus 0.57 19.24 16.80 14.54 0.8 -95.2 
Earth 1.00 10.96 10.40 5.46 3.84 -63.1 
Mars 1.88 5.83 5.50 6.02 1.36 -76.0 
Jupiter | 4346.5 | 2.52x10° 
Saturn 10774.9 | 1.02x10° 
Uranus | 30681.0 | 3.57x107 
Neptune | 60193.2 | 1.82x107 
Pluto 90472.4 | 1.21x10+ 
































It is obvious from Table 3 above that the result of equation (20) appears near to GTR and observation by Hall 
for the first planet, but there is substantial difference between GTR and observation for other planets particularly 
Venus. In the mean time, average percentage of error from prediction using equation (20) and observation (Hall) 
is 7.01%. The numerical prediction for Jovian planets is negligible; though perhaps they could be observed 
provided there will be more sensitive observation methods in the near future. 

It is perhaps also worth noting here, that if we use the expression of quantization of period [3]: 
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T =22.GM.n'/ v (23) 
where v) =@ ¿C= 144km/ sin accordance with Nottale [1]. Inserting this equation (23) into (20), yields: 


D(t) process =A-Qy =V /2GMn°)t / ty (24) 


precess 


or 


T =27/a(t) 


ae = 4T t GMr Ivt (24a) 
This equation (24) and (24a) imply that the anomalous precession of Lense-Thirring type is also quantized. 
Apparently no such an assertion has been made before in the literature. It would be interesting therefore, to 


verify this assertion for giant planets and exoplanets, but this is beyond the scope of the present article. 


precess 


A plausible test using LAGEOS-type satellites 


In this regard, one of the most obvious methods to observe those tad effects as described in this article is using 
LAGEOS-type satellites, which have already been used to verify Lense-Thirring effect of Earth. What is 
presented here is merely an approximation, neglecting higher order effects [12][16][31]. 

Using equation (8c) we could find the rotational effect to satellite orbiting the Earth. Supposed we want to 
measure the precessional period of the inclined orbit period. Then the best way to measure quadrupole moment 
(J2) effect would be to measure the 2 component of the gravity force (8c): 

g =1/r.0V/09=-3GM.a’J,.sin 9.cos.9/r* (25) 

This component of force will apply a torque to the orbital angular momentum and it should be averaged over 

the orbit. This yields a known equation, which is often used in satellite observation: 
@,/@, =—3a°J,.cosi/2r? (26) 


where 1 is the inclination of the satellite orbit with respect to the equatorial plane, a is Earth radius, r is orbit radius 
of the satellite, œ, is the orbit frequency of the satellite, and @,,is the precession frequency of the orbit plane in 


inertial space. Now using LAGEOS satellite data [31] as presented in Table 4: 


Table 4. LAGEOS satellite parameters 





























Parameter Value Unit 
RLAGEOs 12.265x10° M 
ILAGEOS 109.8 o 
TLAGEOs 13673.4 sec 

Ws 4.595x10+ rad/s 
Jn 1.08x107 








Inserting this data into equation (26) yields a known value: 
@, = 0.33756 1°/ day (27) 


which is near enough to the observed LAGEOS precession = 0.343 °/day. 
Now let suppose we want to get an estimate of the effect of Earth kinetic expansion to LAGEOS 
precession. Assuming a solid sphere, we start with a known equation [34]: 


M = AT. P ypnere T? /3 (28) 
where P prere 18 the average density of the ‘equivalent’ solid sphere. For Earth data (Table 1), we get prere 


=5.50x10° gr/m?°. Using the same method with equation (8f), equation (28) could be rewritten as: 
M + AM/At = 47. Py here (r + Ar/ At)’ /3 (29) 


or 





Ar/ At =3/(M + AM 1At) 3 KATP pner) T (30) 


From equation (30) we get dr/dt=13.36 mm/year for Earth. Inserting this value (r+dr/dt) to compute back equation 
(26) yields: 
Ao, =@ @,, =1.41x10°°°/ day =2.558arcsec/ year (31) 


p.nt+l = p,n 
Therefore, provided the aforementioned propositions correspond to the facts, it could be expected to find 
a tad extra precession of LAGEOS -satellite around 2.558 arcsecond/year. To this author’s knowledge this 
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tad effect has not been presented before elsewhere. And also thus far there is no coherent explanation of 
those aforementioned phenomena altogether, except perhaps in [21] and [30]. 
As an alternative to this method, it could be expected to observe Earth gravitational acceleration change 
due to its radius increment. By using equation (28): 
FE) = GM Ir? =42G. Py peel (3 (32) 


From this equation, supposing there is linear radius increment, then we get an expression of the rate of 
change of the gravitational acceleration: 
F(t) = AF / At = 40.G.Py pore (1 + Ar/ At) 1/3- F(t) (33) It would be interesting to find 


observation data to verify or refute this equation. 


Constraint of varying & 


In recent years there is suggestion that unification of the fundamental interactions requires cosmological 
solutions in which low-energy limits of fundamental physical constants vary with time, including o [40][41][42]. 
This assertion began with Dirac’s remark: “the constancy of the fundamental physical constants should be 
checked in an experiment” [42, p.439]. While this has not been widely accepted yet, a plausible way to verify this 
proposition is using celestial quantization method: “offers a possibility to check the variability of the constants by 
studying, for example, lunar and Earth’s secular accelerations, which has been done using satellite data, tidal 
records, and ancient eclipses” [42, p. 441]. 

In this regard, instead of using Nottale’s celestial quantization method [1], alternatively we use Rubcic & 
Rubcic’s assertion [2] that the celestial quantization equation could be related to Planck constant and Planck mass, 


by introducing: 

H'/M = fA (34) 

r, = (fA)’.Mn?/G (35) 
where 

A=h Kam) (36) 
therefore 

r, =(fhlam, )’ Mn? IG (37) 


where f, ñ, @ , mp=2.177x10° kg represents a specific ratio for given system [2], Planck constant, fine structure 

constant (~1/137), and Planck mass, respectively. This alternative expression is quite interesting, particularly if 

compared to Winterberg’s argument of superfluid Planckian phonon-roton as the basic entity of Nature [8]. 
From equation (34) and (36) we get: 


a= fhM (H'm,) (38) 
Using the same method with equation (8f) we get: 

at+a=(M+M).fhi(H'm, ) (39) 
Now dividing both sides of equation (39) with M and @, we get: 

1+a/a=(1+M/M).[fhM (Hm, a] (40) 


From equation (34) and (36) we know that components in the square bracket of the right side of equation (40) 
equals to unity, therefore we conclude by using equation (14): 

àlax MIM =(Mg/ts)/Mg =1/t (41) 

From this viewpoint, we argue that @ varies corresponding to inverse of the epoch of the system in question. 
Supposing the epoch of the Universe is 1.09 Tyr (larger than epoch of the solar system in the previous section, 
2x10'°year), and then from equation (41) we get an estimate of varying @: 

ala ~9.2x10 "year! (42) 

For comparison, other values for varying @ as proposed in the literature are presented in Table 5. 
Alternatively, we could use Sidharth’s original assertion (9) and (14), and by using an equation described in [43]: 

G/G =t; =78.2xdla (43) 
and supposing epoch of the Universe =13.9 Gyr [40d], then we get estimate of @/a =9.2x10" year!, which is 
near to Bahcall et al‘s prediction [40b]. Of course, this subject of varying a@ is not conclusive yet, partly because 
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of large uncertainty in determining the epoch of the Universe, but at least equation (42) could be viewed as an 
alternative constraint based on the celestial quantization method. However, it seems that this proposition of 
varying a from Rubcic & Rubcic’s celestial quantization approach [2] is quite conceivable, particularly from the 
viewpoint of recent suggestion of the invariance and possible time variation of the Planck mass [2b][2c]. 


Table 5. Range of values of varying @ 






































Ref. ala Unit 

Prestage et al. [40, p.31] 3.7x10°4 year! 
Moffat [40d, p.5] 3.8x107° year! 
Ivanchik et al. [42, p.439] 1.9x10!4 year! 
Slyakhter [42a, p.2] 102° year! 
Bahcall et al. [40b,p.520] 2x10! year! 
Bahcall et al. [40b,p.530] 6x10" year! 
Nguyen [40e,p.1] 1.2x10!8 year! 





Discussion 


Allow us to make some more remarks: one of the most relevant being the inconstancy of the "constants". 
Following Dirac's suggestion, the various constants should be validated at period intervals, with no 
assumption that the values will remain exactly the same as they were at the last measurement. 


The fine structure "constant" a , is the most vulnerable to changes in other physical parameters, because it 
is relational. Instrumented observations of the fine structure constant have recorded changes in the value 
starting during the 1980s, and with divergences becoming larger than the original value at the time of its 
inception by Sommerfeld, as time passes. 


The fine structure "constant" is a variable, G is a variable, c is a variable, and if e is variable, the normal 
understandings of the standard physics will be soon be demolished. 


Variations in a are directly related to quantized red-shift, as first pointed to by astrophysicist Halton Arp. 
He also points at quantized changes in gravitation and mass, as related to quasars and active galaxies. It has 
been verified by astrophysical observations, that all galaxies and quasars exhibit quantized red shift. 


But what causes quantized red shift? 


At the core of every galaxy and quasar, there lives an enormous and very complex plasmoid, which 
exhibits periodic episodes of powerful radiations spanning the entire of the E/M spectrum. (There are no 
black holes, anywhere. Those figments are cartoon fantasies produced by Hollywood science to support 
ongoing systematic frauds.) 


E/M sources are also aether sources, and act superluminally. The explosive events of galactic core 
plasmoids produce expanding shells of aether, with abnormal aether density. The various constants, and the 
laws of physics which rely on those "constants" are changed as each aether-density shell expands from the 
center of the galaxy, to the outer edges, changing the local "constants" along the way. 


Where we live, the laws of physics are changing in a gradient manner, along with many of the "constants". 
This gradient of changing values will continue until the entire of the aether shell has passed through our 
solar system. At that point, the constants will remain stable in value, and the physics will change in a 
reliable manner, until the next episode of galactic core-plasma aether ejection occurs and passes through 
where we live. 


The the entire process of changing "constants" will start again, and continue until the shell has passed 
through our location, when the values will once more stabilize and the local physics will again be reliable. 
Until the next galactic core-plasmoid event. 


This is based on current observations of variations in the constants and on quantized red shift, seen in all 
galaxies, as correlated with the SQ originations of all things physical. We are reaching a crescendo of 
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change, including changes of Consciousness and the appearance, ab initio, of new life forms, specially 
constructed to take full advantage of the new conditions which will comprise their environment. (The 
"morphogenic field" of Sheldrake is relevant here.) 


Concluding note 


If physical theories could be regarded as a continuing search to find systematic methods to reduce the 
entropy required to do calculations to minimum; then the fewer free parameters, the better is the method. 
Accordingly, it is shown in this article that some twelve phenomena can be explained using only few free 
parameters, including: 


The Moon is receding from the Earth [20]; 

Earth’s angular velocity decrease (K. Thorne) [19]; 
Planets are receding from the Sun [24]; 

Lense-Thirring effect for inner planets, corresponding to Hall/Newcomb’s observation; 
Celestial orbit prediction in solar system [1][2][3]; 
Exoplanets orbit prediction [1][3]; 

Pioneer-type anomalous acceleration [21]; 

A plausible origin of volcanoes eruption; 

A plausible origin of continental drift effect [29]; 

A plausible origin of spiral motion in spiral nebulae [22]; 
Prediction of extra precession of LAGEOS satellite [31]; 
Prediction of angular velocity decrease of other planets. 





OCOCOCOOOCOCODO 





a 
As a plausible observation test of the propositions described here, it is recommended to measure the 
following phenomena: 


u Lense-Thirring effect of inner planets, compared to spherical kinetic dynamics prediction derived 
herein; 

Annual extra precession of Earth-orbiting LAGEOS-type satellites; 

Receding planets from the Sun; 

Receding satellites from their planets, similar to receding Moon from the Earth; 

Angular velocity decrease of the planets; 

Angular velocity decrease of the Sun. 


Oococvo 





It appears that some existing spacecrafts are already available to do this kind of observation, for instance 
LAGEOS-type satellites [31]. Further refinement of the method as described here could be expected, including 
using ellipsoidal kinetic dynamics [12] or using analogy with neutron star dynamics [32]. Further extensions to 
cosmological scale could also be expected, for instance using some versions of Cartan-Newton theory [38]; or to 
find refinement in predictions related to varying constants. 

All in all, the present article is not intended to rule out the existing methods in the literature to predict Lense- 
Thirring effect, but instead to argue that perhaps the notion of ‘frame dragging’ in GTR [14][16] could be 
explained in terms of dynamical interpretation, through invoking the spherical kinetic dynamics. In this context, 
the dragging effect is induced by the spinning spherical mass to its nearby celestial objects. 
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